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The dynamical behavior of a superconducting quantum interference device (a rf-SQUID) irra- 
diated by a single mode quantized electromagnetic field is theoretically investigated. Treating the 
SQUID as a flux qubit, we analyze the dynamics of the combined system within the low lying energy 
Hilbert subspace both in the asymmetric and in the symmetric SQUID potential configurations. We 
show that the temporal evolution of the system is dominated by an oscillatory behavior character- 
' ized by more than one, generally speaking, incommensurable Rabi frequencies whose expressions 

I are explicitly given. We find that the external parameters may fixed in such a way to realize a 

, control on the dynamical replay of the total system which, for instance, may be forced to exhibit a 

periodic evolution accompanied by the occurrence of an oscillatory disappearance of entanglement 
^ ' between the two subsystems. We demonstrate the possibility of generating quantum maximally 

^ , entangled superpositions of the two macroscopically distinguishable states describing clockwise and 

' counterclockwise supercurrents in the loop. The experimental feasibility of our proposal is briefiy 
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I. INTRODUCTION 



T-H It is of great interest to understand and study the physics of Josephson junction-based devices both for testing 

fundamental properties of quantum mechanics, such as the superposition principle or the occurrence of entangled 
states 0] , and for technological applications in the context of quantum information theory and quantum computing 

Jh ■ In the last decade, rapid developments in the realm of nanotechnologies have made it possible to perform a number 
Ch ' of beautiful and sophisticated experiments at low temperature Q bringing to light the existence in these atom-like 
circuits of many macroscopic quantum phenomena like energy level quantization |Q , macroscopic quantum tunneling 
^ (MQT) and quantum superposition of states ^. More recently, the usefulness of investigating these sohd state 
Q , devices in the context of quantum communication and information theory has been fully recognized. 
O Superconducting Josephson devices may, in fact, be thought of as two-state systems realizing the elementary unit 

^ ^ of quantum information, known as quantum bit or qubit [Q. Moreover, Josephson devices can be scaled up to a 
• ■ large number of qubits and their dynamics may be controlled by externally applied voltages and magnetic fluxes. 
K% ' Superconducting devices like Cooper pair boxes, Josephson junctions or SQUIDs have been thus proposed and used 
^ , as basic elements for the practical realization of quantum gates and chips The relevant macroscopic degree of 
■ - - ' freedom, allowing to store and manipulate quantum information, may be the charge on the island of a Cooper pair 
box or the phase differences at the junction. In the first case the charging energy Ec overcomes the Josephson energy 
Ej. Otherwise, in the opposite regime, the Josephson energy overcomes Ec ||9|. It has been already experimentally 
demonstrated that Cooper pair boxes behave as two level systems which can be coherently controlled 11, 1^ and 
now great efforts are devoted to prove that the same can be done with flux and phase qubits |6|, p^. However, 
successful realization of quantum algorithms critically depends on the ability to entangle quantum states of qubits. 
The optimum would be the realization of a tunable coupling bus. Several coupling mechanisms are possible but the 
natural way of couplin g tw o or more superconducting qubits is through an intermediate resonant LC circuity playing 
the role of a data bus [ p^ Such a resonant LC circuit, describable as a quantum harmonic oscillator, may be in 
principle replaced by the electromagnetic single-mode of a high-Q cavity or by a large-area current biased Josephson 
junction. In all these cases the situation is similar to cavity QED [|l^ (where cavity and atoms play the roles of the 
LC circuit and qubits, respectively) and to ion-trap proposals 
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FIG. 1: Schematics of the superconducting circuits with a flux qubit inductively coupled to a LC resonator modelling a single-mode 
quantized electromagnetic field of a resonant high-Q cavity. 



It is then evidently of interest to study the interaction between a two level solid state system, like an rf-SQUID, 
and an external quantum system like another qubit, a tank circuit or a monochromatic radiation source Jl5| , [T^ , [Tg| , 
pO| , [21] , p2| |2^ . The aim of such investigations is to bring to light the occurrence of entangled states and to 
construct coupling schemes by which the coherent dynamics of the system may be controlled and/or manipulated 
11,11,13. 

In this paper our main scope is to study the dynamics of a flux qubit (an rf-SQUID) coupled to a single mode 
quantized electromagnetic field of a resonant cavity. Confining ourselves to the low-lying energy Hilbert space, we 
prove that the time evolution of the combined system is characterized by the occurrence of entanglement which may 
be controlled in terms of the strength of the coupling and the circuit parameters. In addition we show that the 
dynamics is dominated by an oscillatory behavior traceable back to the existence of a finite set of characteristic Rabi 
frequencies whose expression may be explicitly given. 

The importance of conceiving experimental schemes for realizing quantum superpositions of macroscopically distin- 
guishable states has been quite recently emphasized [ p8| . The main result of this paper is that, appropriately acting 
upon some control parameters, it is possible to guide the rf-SQUID toward coherent maximally entangled combinations 
of two states describing clockwise and counter clockwise supercurrents and then macroscopically distinguishable. 

In section II we describe the physical system under study. Its dynamics in a reduced low-lying energy Hilbert 
space is studied in section III where the main results of this paper are reported. In the last section, we discuss our 
results and we conclude with some remarks about the feasibility of an experiment aimed at verifying our theory in 
the laboratory. 

II. THE QUANTUM CIRCUIT 

In this section we describe in detail the physical system, namely a rf-SQUID coupled to a monochromatic field of a 
high-Q resonant cavity and its hamiltonian model. In figure 1 the electromagnetic single-mode cavity is represented 
as an LC resonator. 

In subsection II A we describe the physical conditions that allow us to consider a SQUID as a two level system, 
that is as a flux qubit. Then we give the hamiltonian model for the cavity field in terms of the combined system 
characteristics and finally, in subsection II C, we consider the inductive coupling between these two subsystems. 

A. The rf-SQUID as material two level system 

Let us begin by describing as usual ||] the rf-SQUID, a superconducting loop interrupted by a Josephson junction 
(figure 2a), as a fictitious particle of mass C and generaHzed coordinate (the magnetic flux in the loop) subjected 
to the washboard potential 

Ui^) = cos (2.^) + (1) 

where C ~ 10~^^ 10"^'^ F is the junction capacitance, L ^ 10 ^ 100 pH the self-inductance of the loop, (j)x an 
externally applied dc flux and 0o = 55 the flux quantum. The Josephson coupling energy Ej is related to the critical 
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FIG. 2: (a) Schematics of a rf-SQUID and (b) of a double rf-SQUID obtained replacing the single Josephson junction with a dc-SQUID, 
namely a superconducting loop interrupted by two Josephson junction. 



supercurrent Ic by Ej — /c(/!'o/27r. 

Taking into account both the kinetic and potential energy, it is then immediate to write the Hamiltonian of the 
system as follows 



where the charge on the junction capacitance Q — —ihd/dcf) and the flux in the loop are canonically conjugate 
operators satisfying the commutation rule [4',Q] = ih. Under the condition [3^ — "^-ttLIc / (j)o > 1 the rf-SQUID is 
hysteretic and the potential U{(j)) may have one or several relative minima. The height of barrier between minima and 
the number of minima depend on the parameter [3l- The form of the potential U{(j)) instead can be tuned by changing 
the external dc magnetic flux (px applied to the loop to the case where the two lowest energy wells are degenerate. 
This case occurs for (p^ — These two degenerate minima correspond to the clockwise and counterclockwise 

sense of rotation of the supercurrent in the loop. The two relative ground states are localized flux states, hereafter 
denoted \L) and \R). For ip^ not exactly coincident with 00/2, C/(0) defines an asymmetric double-well potential 
near (p = (Pq/2 with barrier Vt between these two lowest energy minima. At sufficiently low temperatures, transitions 
between neighboring flux states are dominated by macroscopic quantum tunnelling. However, for high barrier 14, 
tunnelling does not mix the two lowest flux states with the excited states in the two wells. Thus, in the parameter 
regime VJ, ^ fiiuo ^ fc^T ( Huiq being the separation of the first excited state from the ground state in both wells), 
the rf-SQUID effectively behaves as a two state system with reduced Hamiltonian expressible in terms of the Pauli 
matrices and az as follows 



where the basis coincides with the two localized flux states \L) and \R). Here he{(px) = he — 2Ic\/Q{(}l — — 0o/2) 
is the asymmetry of the double-well and A is the tunnelling frequency between the wells. This tunnelling frequency 
can be tuned by changing the height Vb of the barrier which depends on the Josephson coupling energy Ej. Then 
if we replace the junction by a hysteretic dc-SQUID (Fig. 2b), behaving as a J J with tunable critical current, 
Ej may be manipulated by a separately control dc flux (pc- In other words, this modified device, known as the 
double rf-SQUID, behaves as a normal rf-SQUID whose Josephson energy Ej is related to the control flux cpc by 
Ej{(pc) — ^Ico <^os(7r^), Icq being the critical current of each of the two JJs of the dc-SQUID. 

Since in our scheme we need to control A, in what follows, we will investigate a double rf-SQUID with two external 
control fluxes (p^ and (pc when it is exposed to a monochromatic quantized electromagnetic field. 

The hamiltonian of the double rf-SQUID, given also in this case by eq. (H), can be easily cast in diagonal form 




(2) 




(3) 




(4) 



in the basis formed by its eigenstates 



\-)^C^[\R) + 



e + Ve^ + A2 
A 



(5) 
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and 



|+) = C+[|i?) + ^1_!_|L)]. (6) 
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Here = [1 + (^ "^^ j ] ^ normalizing factors of |— ) and |+) respectively and h-sj + is the energy 

difference between their corresponding eigenvalues i?_ = ^ f Ve^ + and i?+ = | Ve^ + A^. 

B. The quantized electromagnetic field 

In this section we model the electromagnetic field of the resonant cavity as a LC resonator. In this way our 
treatment and our results may be easily extended to the case of the coupling between the flux qubit and a real LC 
resonator or a large-area current-biased Josephson junction |jl^, This fact is important because, in order 
to verify experimentally our theoretical predictions, it is more practical to couple the qubit with one of these two 
systems rather than with the single-mode of a resonant cavity. We therefore start by considering the Hamiltonian of 
an LpCp resonator with infinite parallel resistance on resonance and frequency iop = Xj^JLpCp 

2Cp ILp 

where (j)p and Qp play the roles of the magnetic flux and charge operators arising from the quantized electromagnetic 
field and satisfying the commutation rule [(f)p,Qp] = ih. 

As in references jl^, ^ it is possible to relate the resonator operators (j)p and Q_f to bosonic annihilation and 
creation operators a and a^. Defining the flux (pp and the conjugate operator Qp as 



lujpCp 



(a-l-at) 



(8) 



and 



J^ia^a^) (9) 



and substituting these analytical expressions in eq. (0), it is immediate to cast the Hamiltonian of the resonator in 
the form of a standard harmonic oscillator free Hamiltonian 



Hp ^ hujp{a^a+^). (10) 

Here the frequency ujp, choosing Cp ~ 10~^^ F and Lp 0.1 nH, belongs in the range of microwaves {ujp w 
10^^ rad ■ s~^). The eigenf unctions of hamiltonian ( p^ are, of course, harmonic oscillator eigenstates \n) (defined by 
a^a|n) — n\n)) with eigenvalues £"„ = hujp{n + ^). 

If the rf-SQUID is effectively coupled with the quantized mode of an electromagnetic high-Q cavity, we may start 
again from Hamiltonian (0). In order to understand the meaning of Cp in this case, we assume that the rf-SQUID is 
located perpendicularly to the magnetic field and within a distance small compared to the radiation wavelength. In 
such conditions, the vector potential A(x) arising from the electromagnetic field of the cavity mode is approximately 
uniform throughout the region of the device and in international units and in the Coulomb gauge (V • A = 0) it 
assumes the form: 



u being the unit polarization vector, eg the vacuum dielectric constant and V the quantization volume of the field 
mode. Thus the expression for the operator (j)p to be inserted in eq. (|^) may be written down as 



^p = j A-dl=^^—^ia + a^) (12) 
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where the capacitive parameter Cp, given by the following expression 

Cf = e^V j( u • dl^ , (13) 

depends on the field frequency, via the quantization volume V ^ and on the SQUID geometry, via the line integral 
which is taken across a closed circuit 7 inside the SQUID loop. Then, also in this case, exploiting eq. ([l^ ) and the 
properties of conjugation between (pp and Qf-, it is easy to cast Hamiltonian (|^ in the form of a standard harmonic 
oscillator free Hamiltonian (|l0|). 



C. The coupled system 

In view of the assumptions made in the previous subsection, the flux qubit and the LC-resonator modelling the 
monochromatic field can be thought of as coupled together inductively (see figure 1) with a contribution to the total 
hamiltonian given by 

2fc 

Hi = —(t)4)F = B{a + a^)[-f.<j, + l^G^] (14) 

where the constant 



B=^J^ ^ (15) 

L\ 2ujfCf ^/7^T^ 

has the same dimension of h and depends on the coupling strength and the system characteristics. Here the flux 
linkage factor k is assumed of the order of 0.01 in accordance with the current experimental values Thus the 

Hamiltonian describing the combined rf-SQUID-field system can be written down as 

H = Hs+Hf + Hi. (16) 

where Hs and HF^ given by eq. (0) and (nOf) are the Hamiltonians describing the two free subsystems. 



III. THE SYSTEM DYNAMICS IN A REDUCED LOW LYING ENERGY 4x4 HILBERT SUBSPACE 

In view of eqs. (|), (|l^) and (^ the Hamiltonian H may be formally interpreted as that of a two-level 
"atom" interacting quantum mechanically with a monochromatic quantized electromagnetic field. We wish now 
to focuss our attention on a physical situation wherein the operating temperature is T « 10 mK and the field is in 
resonance with the transition between the lowest and the first excited state of the rf-SQUID, that \shijjF = ?iv? + A^. 

As a consequence we drastically simplify the difficult problem of diagonalizing H in its infinite dimensional 
Hilbert space by a truncation procedure which confine ourselves to the 4x4 Hilbert space spanned by the states 
{|0— ), |1— ), |0+), |1+)}. In other words, we investigate the dynamics of the system in the low-lying energy subspace 
of the free Hamiltonian iJp = Hs + Hf generated by the four states, \na) (n = 0, 1; cr = — , -I-) which are the common 
eigenstates of Hs and Hf such that Hs\na) — Ea\na) and HF\na) = huJF{n + ^)\na). 

In this reduced basis and adopting the resonant condition luf ^ \/e^ + A^, Hamiltonian ( p^ ) can be cast in the 
form of the following 4x4 non diagonal matrix: 

/ -Be BA 
_ I -Be huF BA 
~ BA huJF Be 

\ BA Be 2hujF 

We emphasize that the counter-rotating terms of H contribute to this truncated Hilbert space with the presence of 
non diagonal matrix elements of Hji. These matrix elements are those connecting states of the basis having a different 
number n -I- cr of total excitations. 

In the following, we wish to study the dynamics of the coupled matter-radiation system initially prepared in different 
physically meaningful conditions both in the asymmetric (e 7^ 0) and symmetric (e = 0) cases. To this end we need 
to find the eigensolutions of Hn. 
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A. The dynamics of the system in the asymmetric case 

Let us consider the time evolution of the combined rf-SQUID-radiation system initiaUy prepared in the state |0+), 

that is the field in its vacuum state |0) and the SQUID in the first excited state |+) = C+[|i?) + which 
is a superposition of the localized flux states \L) and \B). 

It is worth noting that the manipulation of the shape of the potential and the height of the barrier between the two 
lowest energy wells via the control fluxes (f)x and <^c allows to prepare the SQUID in a prefixed initial condition |9| . 
At the same time, exploiting the currently available experimental techniques in CQED, it is possible to prepare the 
field mode in a Fock state with or 1 photon |30j. Eigenvalues \ui) {i = 1,2,3,4) and eigenvectors Xi of Hr in the 
asymmetric case may be exactly evaluated and they are explicitly given in Appendix A. 

The expansion of |0+)f — exp(— i-^^) |0+) in terms of \u2), jus) and |u4) may be cast in the following form: 

Be 

|0+)= r5 D ^ [VnI{P2 - Qi)Q2 \ui) exp(-iAii) - ^^(Fa + Qi)Q2 \u2) exp(-a2i)+ 

[i^l — i^2)WlW2 

- V^{Pi - Q2)Qi \u3) expi-iXst) + VHI(Pi + Q2)Qi \u4) exp{-iX4t)]. (18) 

Thus |0+) is as in linear superposition, with different weights, of the eigenstates \ui) of Hamiltonian 

We are interested in exploring the ability of the system to periodically come back to the initial state |0+) as well 
as to pass thought the other states |1— ), |0— ) and |1+). The structure of the previous equation and those of eqs. 



( [AlC| )-( A13) makes clear that the time evolution of the system from |0+) involves all the four states |0— ), |1— ), 
|0+) and |1+). However, since the states |0+) and |1— ) are almost degenerate in energy, in our physical situation 
transitions between them are more probable than transitions between the initial state |0+) and the states |0— ) and 
|1+). This fact may be clearly appreciated by evaluating the time evolution of the survival probability Pi{t) of the 
state |0+), that is 

Pi{t) = \{0+\0+)t\^ = S^YS^ + 2SiS2 cos %t + 253^4 cos%i+ 

•> n. n. 



where 



2[S,S3 + S2S4] cos ^^\^^'^ t + 2[S2S, + 8,84] cos ^^^^^0 (19) 



{Pi-P2)QiQ2 



Si = ^iQi~P2) (21) 



82^^{Qi + P2) (22) 



2Be 
P2Q 



2Be 



-{Q2~Pi) (23) 



and 

84 



^(Q2 + Pi). (24) 



It is immediate to construct explicit expressions also for the transition probabilities P2{t) = |(1 — |0+)tp, P3(<) = 
1(0 - |0+)tp and P4{t) = |(1 + |0+)tp to the states |1-), |0-) and |1+) respectively Their explicit analytical 
expressions are given in appendix B. Here we wish to underline that these transition probabilities have the same 
mathematical structure of eq. ([l9|). This means that P2, P3 and P4 like Pi are given by the sum of a constant and 
4 trigonometric time dependent terms with different weights and frequencies. Since these four frequencies uji = 

LU2 = ^f, 1^3 = '^^2h^ ^4 = '^^th^ appearing in the expressions of Pi(i), P2^), Psit) and Piit) are in general 
incommensurable, we find a quasi-periodic behavior wherein a complete exact inversion of the populations between 
the degenerate states 10+) and |1— ) never occurs. The time evolutions of Pi, P2, P3 and P4 are plotted in figure 3. 
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FIG. 3: (a) Survival probability Pi of the state |0+) (dashed line) and transition probabilities P2, (b) P3 and P4 (dashed line) to the 
states |1— ), |0— ) and |1+) respectively, for a system with e 3 • 10-*^" rad ■ , B ^ 2,7 ■ 1Q~^^J ■ s and LOp x lO^^rad ■ s~^. 



A similar quasi-periodic behavior characterizes the time evolution of the system initially prepared in the state 
|0i?) = |0) (E) \R), that is the field in the vacuum state |0) and the rf-SQUID in the localized flux state \R) characterized 
by a well d efine d sens e of c irculation of the supercurrent in the loop. In view of the definitions of |— ) and |+) as well 
as of eqs. ( A14 ) and ( AI6 ), the survival probability of the state |0i?) can be cast in the following form 



Poji{t) = |(0i?|0i?)tP ^W^iy^Wf + 2W1W2 cos + 21^31^4 cos^tH 

— ' ■' n n 



where now 



i=i 



2[WiW3 + W2Wi] cos [Rl_^At + 2[W2Ws. + WiWi] cos l^i±^t} (25) 



W 



2h 



B 



[Pi - P2)QlQ2 



2h 



(26) 



Wi = Q2[eS+{P2 - Qi) + P2A(5_](- 



P2 + Q1 
2BA 



(27) 



W2 = -Q2[eS+{P2 + Qi) + P2AS^]{ 



-P2 



2BA 



Pi 
2Be 



(28) 



and 



W3 = Qi[e6+i-P, + Q2) - PiA,5_](- 



= Qi[e5+{Pi + O2) + PiA(5_](- 



P1+Q2 
2BA 



P1-Q2 
2BA 



5- 



P2 
2Be 



P2 
2Be 



(29) 



(30) 



Here 6± = (\/ + ± e)/2C±-\/e^ + A^. Figure 4 displays this survival probability as well as the transition proba- 
bilities PoLit) = |(0L|0i?)t|2 (dashed line in figure 4a), Pinit) = \{lR\OR)t\'^ and PiL{t) = \{lL\OR)t\^ (dashed line in 
figure 4b) to the states |0L), \1R) and \1L) respectively. Also in this case we underline that the time evolution of all 
these transition probabilities is governed by the 4 characteristic frequencies wi, u}2, W3 and (jJ^ previously defined. 

This leads to a very rich dynamics of the system, characterized by the occurrence of entangled states of the total 
coupled system obtained by the superposition of states with opposite sense of circulation of the supercurrent in the 
loop and a different number of photons in the field (0 or 1). However, due to the fact that these characteristic 
frequencies are not rationally related to each other, it is impossible to restore exactly the initial condition of the 
system. As we will show in the next section, it is possible to find an exact correspondence between these frequencies 
and then a more regular behavior for the total system, in the symmetric case and choosing properly the values of the 
coupling strength k and of the system parameters. 
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FIG. 4: (a) Survival probability Pqh of the state |0-R) and transition probabilities Pq^ (dashed line), (b) Pm and Pu^ (dashed line) to 
the states |0L), \1R) and |1L) respectively, for a system with e ft: 3 ■ 10"'^'' rad ■ B 10~^^J ■ s and iup lO-'^-'^rad • . 

B. The dynamics of the system in the symmetric case: existence of quantum superpositions of clockwise 

and counterclockwise supercurrent states 

In the previous case the asymmetric SQUID potential configuration results from the application of a dc control flux 
(j)^ not exactly equal to 0o/2. In this section we will study the system when (p^ — (/)o/2. This means that the two 
SQUID potential wells have the same height so that e = 0. In such symmetric conditions hamiltonian ( |l7| ) reduces to 
the relatively simpler form 



SA 

hujF BA 

BA hwF 

BA 2nwF 



(31) 



where B must be calculated putting e = in eq. (JTsI). 

Analyzing the structure of matrix ( ^ ) it is not difficult to convince oneself that there exist two dynamically 
separated subspaces, characterized by the frequencies Vli = j-^ up and VL2 — ^'^''+^^ uop, respectively. The first 
subspace is generated by |0+) and |1— ) and the representation of Hr on it is given by the central 2x2 matrix block. 
Such a structure is responsible of the appearance of entanglement in the time evolution of the combined rf-SQUID-field 
system. The matrix elements of H connecting the states |0— ) and generating the second subspace reflect the 
contribution of Counter- Rotating terms in the truncated version of H ^J] . 

The eigenstates of matrix (^l]) assume the simple form 



\Uls 



|0+)] 



(32) 



\U2s 



|0- 



(33) 



\U3s 



B 



|o- 



(34) 



\uis) 



with eigenvalues given by 



(h-VB^ + h^) 



B 



|o-) + |i+)] 



(35) 



Xls = [h- B) LUp 



(36) 



X2s ^{h+B) UUF 



(37) 
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^3s 



{h - UJF (38) 



\is = {h+ y^B^ + h?) UJF. (39) 

Expressions (|3|)-(|3^) and (^)-(|39|) may be immediately derived from those relative to the asymmetric case given 
in Appendix A in the limit for e = 0. The two eigenstates \uis) and |u2s) describe a maximum entangled condition of 
the rf-SQUID and the monochromatic field states, induced by the inductive coupling between them. Rabi oscillations 
between the degenerate states |0+) and |1— ) dominate the dynamical behavior of the system whose time evolution 
may be written down as 

|0+)i = exp [~i\ist/K) + \u2s) exp [-i\2st/h)] (40) 

if the initial condition |0+) is assumed. Also in this case we are interested in exploring the ability of the system to 
periodically come back to the initial state |0+) as well as to pass through the state |1— ). To this end we plot both 
the survival probability Pi{t) = |(0 + |0+)tp = ^(1 + cos2f2ii) of the initial state (solid line in figure 5a) and the 
probability P2{t) = |(1 — |0+)fp — svc? Q,it to find the system in the state |1— ) after a time t (dashed line in figure 
5a). Eq. (|4^) together with figure 5a provides a clear evidence of the existence of coherent Rabi oscillations with 
frequency fii between the states |0+) and |1— ) corresponding to the emission and absorption of a quantum of energy 
hwF by the rf-SQUID. 

Now let us consider the 2x2 subspace in which the dynamics of the system, with respect to the truncated 
Hamiltonian (^l]), is governed only by the Counter Rotating terms with characteristic frequency Vl2. Preparing the 
system in the state |0— ) and considering its time evolution, we easily get 



0-)t = - ...o , „2\ [ / '^'^P {-i^zst/h) -=\uis) exp {-i\ist/h)]. (41) 




Once more we calculate the survival probability P-i{t) = |(0 — |0— )tp of the ground state |0— ) 

td2 At.2 

and the transition probability P4{t) = |(1 + |0— )tp to the state |1+) 

^4(0-^(;^?^[l-cos2M. (43) 

Analyzing the structure of these two expressions, we deduce that the entanglement between the two interacting 
subsystems leads to an oscillatory behavior as before but now, as shown also in figure 5b, we cannot get a complete 
population inversion between the ground state |0— ) and the higher energy excited state 11+). This is due to the 
fact that, in this reduced Hilbert space, processes involving the exchange of two quanta of energy between the two 
subsystems are unlikely. 

Until now we have considered the system initially prepared in a state belonging to one of the two dynamically 
independent 2x2 subspaces. In the following we wish to consider richer physical situations involving both the two 
subspaces at the same time. 

Considering, in fact, as initial condition the state \0R) = |0) (Xi \R), namely the field in its vacuum state |0) and the 
SQUID with a right-hand current in the loop, the system evolves in accordance with the following expression 

|Oi?)t = a{t)\QR) + (3{t)\QL) + j{t)\lR) + S{t)\lL). (44) 

The time dependent parameters appearing in eq. ( ^ ) are linear combinations of trigonometric functions characterized 
by the incommensurable frequencies f2i and ^2- In this case the time evolution of the system is rather similar to that 
obtained in the asymmetric case. Generally speaking, this fact makes it impossible for the system to exactly restore 
its initial condition. Since the ratio ^2/^1 may be controlled by acting upon the parameter B and then, at least, on 
one of the physical parameters appearing in its expression given by eq. (p^, we may wonder on what the dynamical 
properties of the system become in correspondence to special values of B. It is indeed immediate to see that for 



B = h/y/n^ - 1 



(45) 
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(a) (b) 




time (0.1 ns) time (0.1 ns) 

FIG. 5: (a) Survival probability Pi and transition probability P2 (dashed line) at the state |1— ) for a system initially prepared in the 
state |0+). (b) Survival probability P3 and transition probability P4 (dashed line) at the state |1+) for a system initially prepared in the 
state |0— ). In both cases we set e = 0, B Oi 2,7 ■ 10^^^ J ■ s and uip ft: lO^^rad ■ s~^. 



we get ^2 — Qi, where n > 1 is an arbitrary integer. Under such a controllable condition the dynamics of the 
system is dominated by the occurrence of many interesting features. The four time dependent parameters appearing 
in eq. (|4^) assume the form 

, , exp {—iT\/rfi - 1) , . - 1 . , 

a{T) = (cos T + cos nr + I smnr), (46) 

2 n 



^, , exp (— iTVr?— T) / 

Pir) = ^ -{— cos T + cos flT + I 



■ sin nr), 



(47) 



/ \ 6xp {—iT\J rP- — \) . 1 . , 

7(t) = — I (smrH smnr) (48) 

2 n 



and 



^, , .exp {~iT\J V? — 1) , . 1 . , 

5(r) = -i—^ ^(sinr sinnr), (49) 

where r = Slit. Analyzing the time evolution of these four time dependent probability amplitudes, it is not difficult 
to convince oneself that the system comes right back to the initial state \0R) after a time ti = ^ = if n is 

even and after a time ii/2 if n is odd. For this reason eq. ( p5| ) expresses the condition for the occurrence of periodic 
behavior in the dynamics of the combined system. 

Other interesting manifestations in the dynamic of the system occur depending on the parity of the ratio between 
and r^i determined by the fixed value of B in accordance to eq. (|45|). 

We find indeed that, if n is even and always starting from the state \0R), at time ti/2 the combined system reaches 
the factorized state wherein the field is still in its vacuum state |0) and the current in the loop reverses its sense of 
circulation, meaning that the state of the rf-SQUID becomes \L). If n is odd, on the contrary, the probability PoL{t) 
of finding the system in the state |0L) is always less than 1 (see for example for n — 3 the dashed line in figure 6a). 

Moreover, for n even, at times ii/4 and 3ti/4 the system is once more describable in terms of factorized states. For 
example for n = 4 these factorized states may be expressed as 



IMhm = ^"P( y^/') [|0) - z|l)] « (50a) 



= £^pL^^^^[|0)+z|1>]«|_), (50b) 

|— ) being the ground state of the rf-SQUID. This coherent evolution is represented in figure 7a where we plot, for 
n — A, the survival probability PoR{t) (dashed line) of \0R) and the transition probabilities PoL(i), Pinit) (dashed 
line) and PiL{t) to the states |0£), |l_R) and \1L) respectively. 

Analyzing eq. ( ^ ) and eqs. (46)-(|9|) with n = 4, we find that, starting with the field in the vacuum state, at the 



same instants wherein the rf-SQUID get disentangled from the field the probability of finding the rf-SQUID in its 




FIG. 6; (a) Survival probability of the state |0-R) and transition probabilities to the states |0L) (dashed line), (b) |l-f?) and |1L) (dashed 
line) for a system initially prepared in the state \0R). Here, n = 3, e = 0, B = h/\/8 and u>p lO^^rad ■ s~^. 

state I +) exactly vanishes. On the contrary, we may dynamically obtain this condition preparing the system at t = 
in the state |li?). It is indeed easy to prove that in this case, always putting n — 4: in eq. ( p5| ) the state of the total 
system at time ti/A and 3ii/4 becomes once more factorizable assuming the following form: 

|V'i(ti/4)) = £^Ph^^l5/2)[|i) _ ,|o)] « |+), (51a) 

lV'i(3ti/4)) = 2^P(z!^^^[|i)+,|o)]«|+), (51b) 

1+) being the first excited state of the SQUID. Also in this case our theory describes coherent oscillations between the 
states \1R) and \1L), that is the inversion of the supercurrent flow in the loop with period ti/2. This coherent behavior 
of the system assuming \ IR) as initial condition, is represented in figure 7b where, in correspondence to n = 4, we plot 
the survival probability PiB.{t) (dashed line) of |li?) and the transition probabilities PiL{t), PoR{t) (dashed line) and 
-Pot(i) to the states \1L), \0R) and |0i) respectively. It has to be stressed that, on the contrary, when n is odd, the 
previously described oscillations between factorized states and entangled states of the total matter-radiation system 
do not occur. This fact may be fully recognized looking at figure 6, where we plot PoR{t) PoL{t), PiB.{t) and PiL{t) 
assuming n = 3 and Por{0) — 1. We note that a time instant in correspondence of which these four probabilities are 
all equal doesn't exist. 

Remembering that |±) — T 1-^)] and that the states \R) and \L) may be legitimately considered as macro- 

scopically distinguishable states of the rf-SQUID, eq. (^) ((|T])) predicts the generation of a maximally entangled 
Schrodinger cat like state in the dynamics of a rf-SQUID exposed to a single mode quantized electromagnetic field 
when the combined system is prepared in the state \QR) {\1L)). 

The fact of being able to build quantum superpositions of two states describing clockwise and counterclockwise 
supercurrents in the loop, confirms the role of such nanodevices as simple physical systems thanks to which it is 
possible to conceive experiments on fundamental aspects of the quantum theory. 



IV. DISCUSSION AND CONCLUSIVE REMARKS 



In this paper we have investigated the coupled dynamics of a rf-SQUID and a single mode quantized electromagnetic 
field in the reduced 4x4 Hilbert space spanned by the low lying energy states of the uncoupled system. The 
correspondent Hamiltonian model includes contributions from both the rotating and counter-rotating terms and this 
fact turns out to be at the origin of a rich dynamical behavior dominated by Rabi oscillations associated to more than 
one frequency. By construction, our theory is based on a Hamiltonian model containing some external parameters. 
Since they may be easily varied, we have addressed the attractive question of the extent at which this circumstance 
provides an effective tool to get a reasonable control of some aspects of the system dynamics. The analysis reported in 
the paper considers, fixing appropriate resonance conditions, two different cases, the asymmetric and the symmetric 
ones. In both cases the dynamical problem is exactly solved in the truncated Hilbert space finding quasi periodic 
behaviors of the initial state survival probability as well as of some physically meaningful transition probabilities of 
experimental interest. Such quasi periodic temporal evolution reflects the existence of quantum coherent oscillations 
occurring at incommensurable Rabi frequencies. An important difference between the two physical situations under 
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FIG. 7: (a) Survival probability of the state \0R) (dashed line) and transition probabilities to the states |0L), (dashed line) and 

\1L) for a system initially prepared in the state |0-R). (b) Survival probability of the state (dashed line) and transition probabilities 

to the states \1L), \0R) (dashed line) and |0L) for a system initially prepared in the state |l-fi). Here, n = 4, e = 0, B = h/%/TE and 



discussion is that in the symmetric case Hji exhibits two invariant 2x2 subspaces whereas in the asymmetric 
case the time evolution of any initial condition of the total system explores the entire 4x4 Hilbert space. This 
different dynamical behavior has direct remarkable consequences. When indeed the external parameter e measuring 
the height difference between the two minima of U (0) does not vanish, the system evolves intrinsically preventing 
the occurrence of disentangling in correspondence to any possible choice of values for the external parameters. Our 
theory predicts completely different and rich results in the symmetric case stemming from the reducibility into 2x2 
blocks of Hj^ in the truncated Hilbert space. We have indeed proved that the external parameters may fixed in such 
a way to realize a control on the dynamical replay of the total system which, for instance, may be forced to exhibit 
a periodic evolution accompanied by the occurrence of an oscillatory disappearance of entanglement between the 
two subsystems. A relevant result of this paper, commented at the end of the previous section, is the generation of 
quantum superposition of the two macroscopic distinguishable \L) and |-R) of the rf-SQUID. 

These results are new and their value may be further appreciated considering that the realization of our theoretical 
scheme is in the grasp of experimentalists. The several types of SQUID necessary for the observability of our predictions 
are easily fabricated exploiting the well defined trilayer Nb/AlOx/Nb technology. Moreover it is possible to prepare 
and control the state of the rf-SQUID via flux pulses and rf pulses. Finally we may readout the macroscopic flux 
state of the qubit using a suitable magnetometer (essentially an hysteretic dc-SQUID detector) whose experimentally 
measured detection efficiency is of the order of 98 % . One of the most crucial problem related to this kind of 
device is the unavoidable presence of decoherence. The quality of coherence for a two-level system can be qualitatively 
described in terms of the coherence time of a superposition of its states. It is generally accepted that for an 
active decoherence compensation mechanism, T^p must be larger than lOHop, top being the duration of an elementary 
operation of the qubit p2| . In our case the Rabi oscillation frequencies fli and Q2 correspond to characteristic times 
ti « 2.410~^°sec and ^2 ~ 610^^^sec (for n = 4). Moreover, as demonstrated by Cosmelli's group, for a system 
cooled at 5inK and effective resistance i? « 4 SMfi, the decoherence time is approximately of the order of 1/is 
p3[ . Thus in our case we can reasonably believe that it is possible to realize a superconducting device satisfying the 
constraint > lO^top- An open problem for the experimental realization of the physical system is represented by the 
coupling between a qubit and a single-mode of a resonant cavity. We must take into account the typical dimensions 
of the SQUID chip with respect to the size of a high-Q superconducting cavity. In order to bypass this nontrivial 
technical problem we may use an experimental arrangement consisting of a chip placed inside a cavity made by two 
open mirrors [p^ or we may think to integrate the Josephson device and a waveguide in the same chip ||35| . Many 
groups are currently working in this field and we think (and hope) that these techniques will be of common use in the 
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next few years. A most immediate solution is represented by the substitution of the resonant cavity by a LC resonator 
or by a large area current-biased Josephson junction. Several works taking into account this substitution and the fact 
that, as discussed in section III, the hamiltonian model for all these three systems may be expressed in terms of eq. 
(|lo|), make it possible to retain that this experiment may be realized with the currently available technologies. 
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APPENDIX A 

In this appendix we give the analytical expressions for the eigenstates and eigenvalues of Hamiltonian ( |l7| ) in the 
asymmetric case. In order to simplify the notation we introduce the following symbols: 



G = JAB^e^ + h^wl (Al) 



Qi = jAB^Luj, + 2hLUF{hujF - G) (A2) 



Q2 = ^iB^ujj, + 2hajFiruJF + G) (A3) 

Pi = hLUF + G (A4) 

P2 ^huF-G (A5) 
The eigenvalues of the Hamiltonian may be written down as follows: 

X.^hujF-^ (A6) 

A2 = fi^F + ^ (A7) 

A3 = fiu;^ - ^ (A8) 

Xi = hujF + ^. (A9) 

The eigenstates \ui), \u2), {u^} and IU4) relative to the eigenvalues Ai, A2, A3 and A4 respectively assume the 
following form: 

!«■> = l°-> + - 4 + 
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where l/^/nl^ with i = 1,2,3,4, are the normahzing factors satisfying {ui\uj) — 5^ . It is useful to expand the states 



|0— ), |1— ), |0+) and |1+) in terms of |u2), \u^) and |u4). Inverting eqs. ( A10 )-( A13 ) we get 



OA 

= 775 7T\7T7^[^iP2Q2\ui) - ^P2Q2\u2)] - ^.PiQM + V^PiQi|m4)] (A14) 

' /nr(52|wi) - \At^Q2|w2)] - v^Qil^s) + \/^Qi|"4)] (A15) 



(A - P2)giQ2 



= 77^ p',n n ^^^^^ " ^1)^2 - V^{P2 + Qi)Q2 ["2) 

(A — ^2)^1^2 



- V^iPi - Q2)Qi 1^3) + + Q2)0i \ui)] (A16) 



^ ' /T^Q2{PiP2 + Pi-P2Qi+'iB^e^)\ui) + 



2(Pi - P2)Qig2 

- Q2{PlP2 + Pi + P2Q1 + ^B\^) \U2) 

- Qi{Pf + P1P2 - P1Q2 + 4B\^) \U3) + ^ Qi{Pf + P1P2 + P1Q2 + 4B2e2) \ui)]. (A17) 



APPENDIX B 



In this section we give the analytical expressions for the transition probabilities P2 (<) , P3 (t) and Pa (t) to the states 
|1— ), |0— ) and |1+) for a sys te m in the as ymm etric configuration and prepared at t = in the state |0+). 
Exploiting eqs. (Jl^), (A14), (A15) and (A17), it is immediate to written down these transition probabilities as: 



P2(t)-|(i-|o+)tp = 52|^j.2^2riT2 cos^t + 2r3r4 COS^f- 

i=i 



where 



and 



2[rir3 + r2T4] cos iRl-^At + 2[T2n + T^n] COS ^.9l±9At} (Bl) 

Ail L £il L 



T,=Q2^^i-^Bh^ + QiPi) (B2) 
T2 = Q2^^i^Bh^ + Q,P,) (B3) 

n = ^?i^2^(-4i?'^' + Q2P2) (B4) 



T4 = Qi^^(-4i?V-Q2P2). (B5) 
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Also, 

Ps{t) = 1(0 - |0+)t|2 = S^{Y,Z] + cos + 2Z3Z4 cos 

+ 2[ZiZ3 + Z2Z4] cos ^^^2h^'^ t + 2 [^2 ^3 + Z1Z4] cos ^^^^^fp^t} (B6) 



where 



and 



Zi=Q2^^^{-P2-Qi) (B7) 



Z2=-Q2^^{-P2+Ql) (B8) 



Zs = Qi^^{-Pi-Q2) (B9) 



Z^ = Qi^^{-Pi+Q2) (BIO) 



P4W = 1(1 + |0+)(p = 52{y X2 + 2X1X2 cos %t + 2X3X4 cos %i+ 

+ 2[XiX3 + X2X4] cos ^^^^^t + 2[X2X3 + X1X4] cos ^^l±Mt} (Bll) 

where 

= Q2{P2 - Qi) (B12) 

X2 = -Q2(P2 + Qi) (B13) 

^3=Qi(-Pi+Q2) (B14) 

X4 = gi(Pi+Q2). (B15) 
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